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O ; ABSTRACT 

' , Time evolution of tunneliner in thermal medium is examined using the real-time 

Ph. 
^ • semiclassical formalism previously developed. Effect of anharmonic terms in the 

potential well is shown to give a new mechanism of resonance enhanced tunneling. 

K> , If the friction from environment is small enough, this mechanism may give a very 

c^ ' large enhancement for the tunneling rate. The case of the asymmetric wine bottle 

potential is worked out in detail. 



In our previous paper |1|] we formulated within the semiclassical framework how 
to deal with the real-time dynamics of tunneling that takes place in thermal medium. 
In the present work we further analyze the problem and discover a new mechanism 
of enhanced tunneling caused by anharmonic terms in the potential well. The result 
may have great relevance to the old (once failed) scenario of inflation based on the 
first order phase transition [^, and to the scenario of electroweak baryogenesis |0|. 

We start from recapitulating our previous result [jl|. The real-time dynamics 
of tunneling in thermal medium is studied in the standard model [Q,]^, of en- 
vironment. Its interaction with a one dimensional system which we describe by a 
potential V{q) is of a bilinear form, q J du c{uj)Q{uj) . Here Q{uj) is the environ- 
ment oscillator coordinate of frequency u and c{uj) gives a coupling strength of the 
system-environment interaction. With the total system thus specified, dynamics is 
given by the quantum equation of motion, 

+ -r = - diu c{iu)Q{uj) , X, +^ Qiuj) = -c{uj)q. (1) 



dt'^ dq Juic ' dt"^ 

Quantum Langevin equation is derived ^ by eliminating the environment vari- 
able Q{u! ,t); 

d'^q dV /■* / ,,,„,, 

where Fqlt) is linear in initial environment values, Qi{uj) and Pi{uj), 

roo ( P(uj) \ 

FQ{t) = - / dujc{uj) Qi{uj) cos{Lut) + ^^^ sm{ijt) , (3) 

and obeys the correlation formula, 

({FQ(r),FQ(s)}+)env = / duj r{uj) COS uj (t - s) coth{—-) = a r{t- s), (4) 

with r{uj) = c'^{uj)/{2uj) and /? = 1/T the inverse temperature. The kernel function 
aj in eq.(l) is given by ai{t) = — J^ dujr{ijj) sin(u;t) . The combination, auit) + 
iaj{t), is a sum of the real-time thermal Green's function, added with the weight 

An often used simplification is the local, Ohmic approximation taking r{uj) = 
rju/n with Uc = 0, which amounts to a/(r) = 6uj^6{t) +?7(5'(r) . This gives the local 
version of Langevin equation, 

d'^q dV ^ 2 dq ^ , , 



The parameter duj"^ is interpreted as a potential renormalization or a mass renormal- 
ization in the field theory analogy, since by changing the bare frequency parameter to 
the renormalized uj\ the term duo"^ q is cancelled by the counter term in the potential. 
On the other hand, rj is the Ohmic friction coefficient. This local approximation 
breaks down both at early and at late times |p, but it is useful in many other cases. 
The crucial observation |I| is that in the semiclassical approximation the Wigner 
function has the solution in terms of classical solutions; 

fw{q,P, Q ,P;t) = J dqidpi J dQidPi fw{qi,Pi , Qi , Pi) 

■S {q - gel) ^ [p - Pel) 5 [Q - Qd) 5 (P - Pd) . (6) 

Here qc\{qi ,Pi ,Qi ,Pi;t) etc. are the solution of (|1|), taken as the classical equation 
with the specified initial condition. We take in the present work an uncorrelated 
initial state of the form p^^^ = p^*) Pn ■ The reduced Wigner function is defined 
after the environment variable {Q{uj) ,P{uj)) integration; 



fw\(l^P'^t) = J dqidpi fi^^giqi , Pi) K{q ,p ,qi ,pi;t) , (7) 

K{q ,p,qi,Pi;t) = J dQidPi fw ,Q{Qi , Pi) Hq - Qd) Hp - Pa) ■ (8) 

The next important point is expansion of the classical solution g^i in the environ- 
ment initial variables {Qi , Pj); 

gci(t) ^ qS\t) + j du[ QMqf\u , t) + P,{u) qiPiu; , t) } . (9) 

The truncation up to the first order term in {Qi{uj) , Pi{uj)) should be adequate at low 
temperatures. With the help of the familiar Fourier formula for the delta function in 
eq.(|D, one derives, using the Gaussian thermal function for the initial fw ,QiQi ; -Pj); 
an integral transform of the Wigner function, fw^q{qi,Pi) — > fw ■ The resulting 
kernel function is given by [|I| 



K{q ,p,qi,Pi;t) = — exp 

ZTT 



^ \ P-Pcl 



(10) 



The matrix elements of {J' ^)ij = lij are given as follows. First, 

1 f^ ^ ^^PCO 1 ,_,^ ^,,2 



hi = - duJCoi\i^-\z{uj,t)Y, (11) 

and I22 is given by a similar integral, replacing z{u! ,t) in In by z(u!,t), while 
/12 = 2^- Here z{uj ,t) = q^^i^u ,t) + iujq^^i\u ,t) , and z{u ,t) = p[.P{uj ,t) + 
iujpiP{uj,t). 



Quantities that appear in the integral transform are determined by solving differ- 
ential equations; the homogeneous Langevin equation for g^j and an inhomogeneous 
linear equation for z{u! ,t) and z{u! ,t), 



\dq ) (0) Jo 



(12) 



d''z{u , t) ^ fdH^ I z{uj,t) + 2 I ds aj{t - s)z{iu , s) = - c(cu)e*"* . (13) 



dt^ \dq^ 



9, 



(0) 
cl 



A similar equation as for z{uj ,t) holds for z{uj ,t) . The initial condition is 
q^\t = 0) = q,, p^\t = 0)=p,, ziu,t = 0) = 0, i(u;,t = 0) = 0. 

The physical picture underlying the formula for the integral transform, eq.(0) 
along with (|10D , should be evident; the probability at a phase space point {q,p) is 
dominated by the semiclassical trajectory q^^ (environment effect of dissipation being 
included for its determination by eq.(|T^)) reaching (g ,p) from an initial point (g^ ,Pi) 
whose contribution is weighed by the quantum mechanical probability f^/ initially 
given. The contributing trajectory is broadened by the environment interaction with 
the width factor Jlij- The quantity In given by ([TT| ) and (p!3D, for instance, is equal 



to (g — gel )^ 5 ^^ environment driven fluctuation under the stochastic force FQ{t). 
The tunneling potential is divided into two regions separated at the position 
q = qB oi the barrier top. In the present work we shall assume that the potential 
is very steep at both ends; ^(±oo) = oo . The tunneling rate, from the inner 
region at g < g^ called here the potential well into the outer region at g > g^, is 
an important measure of tunneling phenomena and is given by the flux at g = g^; 
P{t) = — I{qB , t) , which is equal to 

- / dq^dp^ /« (g. , p^) (^J) + ^ (^B - ,S' ) ) ^ exp [ - ^^-;f^' 

The width factor ^/Tu thus determines the contributing region of g^j by how much 
away it is from g^. On the other hand, the tunneling probability into the overbarrier 
region at g > x is given by 



dqidpi f^' {qi , Pi) / du ^— ^^ exp 



{u - qf Y 



2/ 
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(15) 



In both of the quantities, P{t) and P(g_B ,^) the tunneling probability into g > g^, 
it is essential to estimate how g^i and In varies with time. 



We shall assume an initial state localized in the potential well so that the domi- 
nant contribution in the {qi ,Pi) phase space integration is restricted to g^ < qs- Let 
us first consider the harmonic approximation for the potential well region given by 
V{q) + ^ 6uj^ g^ ~ I cUq q"^ near the bottom of the well at g = 0, and use the Ohmic 
friction. The approximate solution is given by 



g ^^ = coscDot + -^ sinupt e""*/^ q, + -^^ e-^"' p, , (16) 

V ZUOq ) Ijjp 

c(c<j) 
ijj^ — {jJq— tujr] 

V 2c^o 2cJo / 



„2 



using Up = \Juop — ^ . We assume a small friction, rj <^ uop . 
Near uj = Cop this function is approximately 

z(oo,t)^ '""^"^^ , fte^^"*- — sincZ^oO e~^'/\ (18) 

UJ + ujp — iri/2 \ Up J 

This formula is valid at t < 0(1/?]]. The appearance of the linear t term is a resonance 
effect. The resonance roughly contributes to /ii(t) by the amount, rjt^ e~^* x a 
smooth u integral which is cutoff by a physical frequency scale. Thus, the width 
factor 111 initially increases with time until the time scale of order l/rj. 

At large times the classical trajectory g^j asymptotically approaches towards 
g = 0, the local minimum of the potential. The width factor behaves as 

Jn(t)=/ii(oo)+0[e-''*/2], (19) 

In oo ^ — + - ,^ , + tt4 ^ • 20 

2ujp ojp e'^o/i _ 1 3 ^g 

The asymptotic value of /ii(oo) has the familiar zero point fiuctuation of harmonic 
oscillator and in the last term the dominant finite temperature correction, valid for 
this Ohmic model at T <^ ujp. In any event the probability rate P{t) finally decreases 
to zero, along with 

p^r+^fe-.i^-o. (21) 

Moreover, the final tunneling probability P{qB , oo) has a finite value, and typically 
is very small for a large potential barrier. For instance, for the asymmetric wine 
bottle potential later discussed. 



^"-~' 4 1/2^^" '"■'""• P^' 



with Vh the barrier height much smaller than cjq. This poses a curious question; it 
appears that decay of a prepared metastable state localized in the potential well is 
never completed. 

This simple picture is however valid only when one ignores anharmonic terms in 
the tunneling potential, but they must be there in order to give any realistic tunneling 
potential. The most important point for the present work is effect of anharmonic 
terms in the equation for the fluctuation z{uj ,t). Presence of anharmonic terms 
gives a non-trivial periodicity in the coefficient function (^-r) (o) for ([13D , assuming 
a small friction r] <^ uoq. For a very small friction the anharmonic term becomes 
important both for completion of the decay and the new mechanism of resonance 
enhanced tunneling. In the rest of this paper we shall discuss the mechanism of 
resonance enhancement. 

One might naively expect that the homogeneous part of the z{(jj ,t) solution 
exhibits the well known parametric resonance 0, if the relevant parameter in the 
periodic coefficient function (^-r) (o) feUs in the instability band. But this is not 

"cl 

what happens, as will be shown. Thus, unbounded exponential growth of -\//n does 
not take place. On the other hand, the power-law growth is observed, as clearly seen 
in numerical computation. Moreover, we find that our z{uj ,t) solution belongs to 
the boundary between stability and instability bands. At the resonance frequency 
the enhancement factor due to the boundary effect is much larger than what one 



might expect from the harmonic case (pISl). We shall interpret this phenomenon 



as influenced in a subtle way by the parametric resonance, although it is not the 
parametric resonance itself. 

Observe first that if q^ {t] C[i,Pi) is a solution for the classical, homogeneous 
Langevin equation, then q^i {t + e; qi,Pi) for an arbitrary constant time shift e 
is also a solution. Hence the time derivative qli {t; qi,Pi) = Pc\ i^', qi,Pi) is a 
homogeneous solution for the linearized z{uj ,t) equation. Using this homogeneous 
solution, one can readily find the inhomogeneous solution for z{uj ,t) with the initial 
condition, z{iu , 0) = z{u , 0) = . Corresponding to the classical motion q^i of given 
initial condition (g^ ,pi), this solution is 

z{u,t) = -c{u;)pf{t)fdt' {p^\t')y^ f cit"e-^'(*'-*")p(p(t")e*"*". (23) 

It is best to discuss the resonance enhanced tunneling mechanism in concrete 
examples. We take the asymmetric wine bottle potential as illustrated in Fig.l, 



which is described in the well and its vicinity region by 

V{q)^\{q^-2qBqf. (24) 

The curvature parameters at two externa of g = and q = qs are Uq = 2Xq\ , uj\ = 
\q% , and the barrier height seen from the bottom of the well is Vh = \q% = uJ^q^/S . 
The condition for a large barrier, Vh ^ t^o, implies \fXq% ^ 4-\/2 . The classical g^j 
equation in the Ohmic approximation is written using rescaled variables, y = q\.\ /qs 
and T = ujq t/2, 

l/" + 2/(y-l)(y-2) + ^y' = 0, (25) 

where y' = dy/dr. 

To give an idea of the magnitude of initial amplitudes, y{0) and y'{0), let us 
assume an initial thermal state trapped in the (hypothetical) harmonic potential well, 
taking the same temperature as the environment. The magnitude of initial values is 



then on the average y qf = y'coth ^/(2u;o) ~ vl/(2^o) , y pf ^ y tUo/2 • This 
gives 

y qf yPi i /^ 



qs ^oqs 

.(0) 



(26) 



The classical solution p^j that appears in the z{iu ,t) solution (p^) is given by 
the Jacobi's elliptic function for the symmetric double well, the situation relevant to 
the case of a small friction. This function has the fundamental period given by 



1 du 



UJo a/I + V^ ^° J(1-m2)(1_A;2u2) ' V + V^ ' 



'^ (27) 



where e = Ei/Vh , Ei = pj/2 + lJq qf/2 . The corresponding fundamental frequency 
u^: is 27t/T. The parameter k, hence e is a measure of the anharmonicity. For 
instance, E{qi ,pi) = Vh/A corresponds to A; ~ 0.80, a value giving cj* different from 
the harmonic frequency uoq at the potential bottom; uo^, ^ 0.95 x ujq. These formulas 
of the period are used to later determine the band structure associated with the 
z{u ,t) equation. 

We numerically integrated the coupled g^i and z(ij ,t) equations. In the to inte- 
gral (0) for the width factor In the largest contribution is found to come from the 
resonance ai u = u^, then contribution from higher harmonics at nu^ follows. Thus, 
phenomenon of a non-linear resonance occurs. In Fig. 2 we show the width factor 
as a function of time for a few choices of the initial data (g^ ,pi). For simplicity an 



equi-partition of the kinetic and the potential energy pf/2 = uJq qf/2 is assumed to 
reduce relevant parameter dependence. The result for the width factor is plotted for 
a few choices of the initial energy. At an intermediate time of ©[l/r^] the width factor 
becomes maximal at its value much larger than in the harmonic case. The decrease 
at late times is due to the friction; we explicitly checked that for the zero friction, 
the fluctuation \z{uj ,t)/c{uj)\'^ increases in time without a bound, with an averaged 
time power close to 4 at the resonance. Note that even if the effect of the friction 
is turned off, there exists an important environment effect here; the environment 
interaction drives the non-linear resonance oscillation. 

The unbounded power-law increase suggests that the relevant parameter in the 
homogeneous z{uj ,t) equation falls in the boundary between the stability and the 
instability bands, since otherwise it either is bounded or grows exponentially. We 
indeed checked this numerically by arbitrarily changing the parameters in the z{uj ,t) 
equation. We used the homogeneous equation introducing two arbitrary parameters 
{h,e); 

z"+(^h-2e{y-^-y')^z = 0. (28) 

Unlike a single, linear equation for the Mathieu type our problem is a coupled, non- 



linear system, with y to be determined by eq. f^Sf ) for rj = 0. We must further 



assume some initial values for y{0) , y'{0). The standard algorithm that determines 



the stability and the instability band is used |T0|. The parameter set {h ,6) = {4 , 3) 
corresponds to our homogeneous z{uj ,t) equation for the zero friction. This case is 
seen right on the boundary line of the two bands as depicted in Fig.3. We checked 
that our system is always on the boundary line for all choices of y{0) , y'{0) we 
computed for. 

A rough analytic understanding of the growth of the width factor seems to be 
as follows. In the explicit formula for z{uj ,t), eq. (P3|) , one may replace the quantity 
[Pel {t')) in the integrand by some constant average value. For a small friction it 
is then easy to obtain a linearly growing z{uj ,t) off the resonance and a quadratically 
growing z{uj , t) right on the resonance as a function of time in the range t < 0[l/i]]. 
This time dependence is what is observed in numerical computation. 

Computation of the tunneling rate P{t) is a demanding task of numerical cal- 
culation, and it is left to future work. But a general trend is expected already in 
the following simplified computation. The most important part for P{t) is two com- 
peting exponential factors, the one for the initial state and the other for the kernel 



factor; 

(9B - g?^ ? 



A = exp 



,,,^^p1±M 



X exp 



2/ 
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(29) 



We plotted in Fig. 4 the time evolution of this product. The time averaged evolution 
(averaged over a short time scale of several times the fundamental period) shows a 
rapid rise of the product factor, reflecting the initial behavior, In oc t^. The product 
then reaches some maximum at a time of order l/rj. The asymptotic time limit 
coincides with that given by the harmonic approximation. As a numerical guide 
note that e-«s/(2-fii) = f,-{2Vjcvo){ujohi/2)-^ ^ g- 20(0.0/11/2)-! ^^ y^ ^ ^Qcuq. 

In Fig. 5 we show for a few choices of the friction the maximum value of A as a 
function of the initial energy Ei, assuming the equi-partition pi = uo^qi. Remarkably, 
the largest contribution comes, not from the dominant initial component near the 
zero point energy, rather from the initially suppressed excited component. The first 
exponent in (pQ]) is in proportion to i?,, while the second one goes roughly like -Ej"^, 
hence a maximum may appear somewhere away from the lowest energy of ii^j = co'o/2. 
In this example of tj/ujo = 0.0025 the maximum product factor is of order 10^'^ at 
Ei ^ A X uJf)/2, 6 orders of magnitudes larger than what one expects from the 
lowest energy state and also the asymptotic value of order 10"'^^. A large value of 
the product factor of order unity suggests an interesting possibility of a rapid and 
violent termination of the tunneling. 

We end with some speculative comments on impact of the resonance enhanced 
tunneling when applied to cosmology. For this part of discussion we assume that the 
non-linear resonance we found in the present work terminates the tunneling associ- 
ated with the asymmetric wine bottle type of potential. The tunneling phenomenon 
is clearly related to termination of the first order phase transition, if there exists a 
prolonged period of cooling providing a chance that the universe, or a part of it, is 
trapped in a supercooled metastable state localized, say near g = in our quantum 
mechanical terminology. The coordinate q here should be understood as the order 
parameter of the phase transition, the homogeneous Higgs field in the case of GUT 
and electroweak 0] phase transition. The environment in this case is made of 
various forms of matter fields, the quark and the lepton fields, and also the gauge 
and the Higgs particle. In the usual weak coupling theory the friction containing 
typically small Yukawa couplings is small compared to the typical energy scale of 
the system, the Higgs mass. Thus the small dissipation seems a good approxima- 
tion. The low temperature approximation also works if T <^ the Higgs mass. After 



many periods of the periodic Higgs motion (the period ~ 27r/(Higgs mass) ) the res- 
onance amphfication is expected to lead to an outburst of tunnehng. This may occur 
stochastically in different parts of the universe, leaving behind nucleated bubbles, al- 
though it is not clear how different bubbles subsequently merge. Since the nucleation 
occurs presumably within a short time in all parts of the universe, it appears that 
the first order phase transition is terminated abruptly. This highly non-equilibrium 
incident should much help the electroweak baryogenesis. 

The other possibility of completing the phase transition is signaled by the change 
of the potential (more properly of the free energy) due to the temperature change so 
that the once local minimum at g = becomes an unstable, local maximum of the 
potential. Under a normal circumstance the rate of the temperature change is given 
by the Hubble parameter H, hence if the Hubble parameter is small enough, like 
H < f] (the friction), completion of the phase transition via the resonance enhanced 
tunneling most likely takes place. This condition seems to be satisfied in the two 
cosmological examples given. 

In summary, we gave using the real-time formalism of tunneling dynamics a new 
mechanism of resonance enhanced tunneling. 
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Figure caption 

Fig.l 

A schematic form of asymmetric wine bottle potential. A case of large barrier 
Vh = 10 ujq is depicted. The barrier top is at qb = 4:^/5 /^/uJo. 



Fig.2 

Time evolution of the width factor In. A time average over uJoAt/2 = 20 of the 
dimensionless quantity ujoIu/2 is plotted as a function of time xujq/2. Examples of 
initial states of energy ii^j = (1 , 2 , 5) x ujq/2 {ujq/2 being the ground state energy 
in the potential well) are compared to the harmonic case. Assumed parameters are 
Vh = lOcJo , ri/ujQ = 0.0025 , T/ujq = 1/20 , fi(cutoff)/c<;o = 25 . 

Fig.3 

Band structure of stability (unshaded) and instability (shaded) for the modified 
fluctuation equation. The defining equation is (|28|), along with (^), with the initial 
condition of y{0) = y'{0) = 0.2 . The point marked by the cross has {h ,9) = (4 , 3) 
corresponding to our system. 

Fig.4 

Time averaged evolution of the product A of two exponential factors for the 
tunneling rate. The same parameter set as in Fig.2 is taken. In the inlet the case of 
Ei = 5 X to'o/2 is depicted in the linear scale. 

Fig.5 

Maximal value of the product A of two exponential factors for the tunneling rate 
is plotted as a function of the initial energy for a few choices of the friction t]. The 
other parameters are taken the same as in Fig.2. 



12 



Fig. 1 
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